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Using the finite Fourier transform, we introduce a generalization of Pauli-spin matrices for d- 
dimensional spaces, and the resulting set of unitary matrices S (d) is a basis for d x d matrices. If 
N = di x d 2 x ■ ■ ■ x d b and H [N] = (g)# [dfcl , we give a sufficient condition for separability of a 
density matrix p relative to the _ff' dfc l in terms of the L\ norm of the spin coefficients of p. Since 
the spin representation depends on the form of the tensor product, the theory applies to both full 
and partial separability on a given space H^- N K It follows from this result that for a prescribed form 
of separability, there is always a neighborhood of the normalized identity in which every density 
matrix is separable. We also show that for every prime p and n > 1 the generalized Werner density 
matrix W^ p ' (s) is fully separable if and only if s < (l +p n_1 ) . 
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r*"** '■ I. INTRODUCTION 

One of the predictions of quantum mechanics is that spatially separated components of a system can be entangled, 
■^j- \ The consequent prediction of non-classical correlations among the separated components of a quantum system has 
led to critiques of the foundations of quantum mechanics, as in the famous Einstein, Podolsky, Rosen paper [Q, 
and to experiments that have confirmed the predicted non-classical correlations, as in Q. Interest in entangled 
systems has been heightened by proposed applications in quantum computation, for example fsj, and in quantum 
communication, as exemplified most dramatically by teleportation Q. As a result there have been many publications 
which have examined various aspects of entanglement, its measurement, and its use in quantum communication such 
as references to mention only a few recent papers. 
| In this paper we shall be interested in the separability properties of quantum systems in states defined on finite 
dimensional Hilbert spaces H = Hi ® ■ ■ ■ <g> H n , where the H^ denote the Hilbcrt spaces of the subsystems. A state 
specified by a density matrix p is said to be completely separable on H if it is a convex combination of tensor products: 

p = Y J P{a)p {1) {a)®---® P ^{a), (1) 
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where p( k > (a) is a density matrix on Hj~. Since the same p can have different convex representations, it has proven 
difficult to determine generally applicable operational conditions for separability, and determining such conditions is 
one of the motivations for this paper. It is also possible to have different types of separability by allowing sets of the 
subsystems to be entangled, cf. P,|fO|| , and one can describe a lattice of levels of separability. The theory we develop 
here applies to all of these various definitions of separability. 

A necessary condition for separability is that the partial transpose p Tr of a state p should be a state |tl]]. If we 
represent p as a matrix, this means that if p = (pj 1 ...j n ,ki...k nt ) then (taking r = 1) 

is also a density matrix. It is easy to confirm that if a density matrix is separable, its partial transposes are also 
separable, but it has been shown [EJJ that the converse is true only in the 2 ® 2 and 2 ® 3 cases. In the proof of 
this last result a necessary and sufficient criterion for separability was established, but there seems to be no 
operational way of using this criterion as a general tool. Other studies of separability, such as those in |p|,|f3|-p^| have 
found operationally useful necessary conditions and sufficient conditions for classes of densities or for special cases, 
but no general sufficient conditions with a breadth of applicability analogous to that of the Peres condition. 
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Broadly speaking, necessary conditions tend to be described in the computational basis while sufficient conditions 
for 2-level systems tend to be described in terms of the Pauli spin basis. That observation motivated the derivation 
of a change-of-basis formula in |l6| which facilitates the strategy of checking whether necessary conditions derived in 
the computational basis are sufficient by using the (real) Pauli spin basis. This approach leads to general sufficient 
conditions for full separability which essentially give the condition in H] as a corollary and also leads to necessary 
and sufficient conditions for full separability of a parametrized family of n-qubit densities which all satisfy the Peres 
condition. The difficulty with extending this approach to d-level systems is that the generally accepted definition of 
spin matrices as generators of rotations does not capture the computationally useful features of the Pauli matrices 
when d > 3. One of the basic purposes of this paper is to propose a general definition of d-level spin matrices which 
possess many of those computational properties. 

The Pauli matrices are special in that they are both Hermitian and unitary, and together with the identity matrix 
<7o they form a basis of the set of 2 x 2 matrices. Our strategy is to generalize the role of the Pauli matrices as 
a basis of unitary matrices at the expense of Hcrmiticity. We show the applicability of these proposed d-level spin 
matrices and in the appendix examine the d — 3 case in some detail, identifying properties analogous to those of the 
Pauli matrices. We also define directly a general characterization of certain classes of trace one projections of d-level 
systems. We use those projections to establish necessary and sufficient conditions for full separability of generalized 
Werner densities composed of any number n of d-dimensional subsystems for any prime d. In addition, we establish 
a general sufficient condition for full or partial separability of densities of any dimension. Analogous results for full 
separability were obtained in |f6| for d = 2 by essentially the same methodology. 

Other authors jlTj have used a different set of operators in the d = 3 case and some separability results were obtained 
recently in []l5|| . Our proposed class is different, and we show that stronger separability results can be obtained using 
these matrices and the strategy developed in E3j. 



II. A NECESSARY CONDITION 



As mentioned above, the Peres partial transpose condition is a general necessary condition for separability pT| . In 
[p~6|| a weaker but useful condition was derived using the Cauchy-Schwarz inequality and has the following application. 
Suppose j = ji . . . j n and k = k\ . . . k n differ in each component: j r k r . Let u and v be indices with u r =/= v r and 
{u ri v r } = {j r , k r } . Then for fully separable states p 

(VPjlVPKk) > \Pu,v\ , (3) 

where p is written as a matrix in the computational basis defined by the tensor products of \ji) (ki \ , 1 < i < n. As 
an application, consider the following generalization of the Werner density matrix fl8[| on the N = d n dimensional 
Hilbert space H 

w m (s) = i_i 7 + S . T 

v ' d n 

where / is the identity and r is the projection defined by the state 

^]) = -L(| ...0) + |l...l)+... + |(d-l)...(d-l))). (4) 

(In the sequel we let k denote the repeated index k . . .k.) In the computational basis W^^ (s) equals (^jsr + g ) 
when j is in jfc : < k < d| and equals otherwise. The only non-zero off-diagonal elements are 

4. Choosing j and k appropriately in (||), we have the necessary condition 1 > s (l + d" -1 ). To show that this 
condition is also sufficient, we will use the spin representation to prove (s) is fully separable when d is prime 

and s = (f + d" _1 ) . In order to do that, however, we first need to define the spin representation. 



III. COMPUTATIONAL AND SPIN BASES 

Let if [ w l denote an A^-dimensional Hilbert space where N — d\ x d 2 x ■ ■ • x db- In this section we define different bases 
for N x N matrices on based on different representations of as a tensor product space, and the discussion 
is purely mathematical. In the applications we will be concerned with a specific representation = <§) a=1 H^ da ^ 
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and with the corresponding separability properties of densities on H^ N \ The bases used will depend on the order of 
the tensor product as will the representation of a density matrix as p = p\ (g> • • • <g> pb, the tensor product of densities 
Pk on the H^ dk \ For example, we might want to examine separability of a density matrix on ffl 90 ! = <g> 
using matrices consistent with that tensor product. In a subsequent application one might want to represent 
the tensor product of spin 1/2 and spin 1 particles, i.e. a tensor product of and W 3 ', and the order of the 
sub-tensor product shouldn't affect the theory. We confirm that assertion by showing that permuting the order of 
a tensor product corresponds to a conjugation operation and thus that the theory is generally applicable with only 
notational changes for particular applications. 

Lemma 1 Let N = d% x d 2 x • ■ ■ x db and suppose M is an N x N matrix with M = ® •• • ® C^ b \ where the 
C (fc) are d k x d k matrices. Let a denote a permutation of {1, ... , b} and let M a = C (<t(1)) ® • • • ® COW). Then for 
all such C^'s there is a permutation matrix Q a such that Qa-M^Q" 1 = M. 

Proof. If M {j, k) = (j\,k\) (Jb, kb), the index j corresponds to the ordered 6-tuple (ji, . . . ,jb) and j 

is uniquely defined by j = ji [d% x • • • x db) + ji {d-& x • • • x db) + • • ■ + jb and similarly for k. Let j CT correspond to 
(io-(i)i ■ ■ ■ >j(r(b)J an d define the permutation matrix Q a by Q a (j, s) = S (j a , s). Then Q^ 1 — Q l a and 

Q„M a Q- x (j, k) = M a {j^K) = k a(1) ) ( Mb) ,K {b) ) = M (j, k) , 

completing the proof. 

The motivation for this work comes from the identification of the 2x2 Hadamard matrix as a key tool in working 
with 2-level systems. Specifically, suppose p = | (er + a m ), where a m = J2 m j a j anc ^ a i = a xi a i ~ u v an< ^ ^3 = °2 
are the usual Pauli matrices. Then the coefficients in the spin basis are related to those in the computational basis 
by the 2x2 Hadamard matrix: 



1 mi \ _ / 1 1 \ f poo Poi 
m 3 ~im 2 J \ 1 -1 J \ pn pio 



(5) 



The matrices in the two bases are connected in a similar fashion: 

(6) 
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-^0,0 ^0,1 

£■1,1 £-1,0 

where the matrices Ej t k — \j)(k\ define the computational basis. Note that a systematic application of these re- 
lationships requires both the use of the real Pauli matrices and a reindexing of both the Pauli matrices and the 
computational basis matrices to conform to the matrix notation. 

The Hadamard matrix is the 2x2 Fourier transform, and we extend the idea of (|^) to d x d matrices. The adjusted 
basis A = {Aj^, < j,k < d} is the set oidxd matrices defined by A,-.fc = Ejj + k, where + denotes addition modulo d, 
and we define the "spin" matrices S = {Sj t k, < j,k < d} using the analogue of (^) and the finite Fourier transform. 
Thus (S) = F ■ (A) where F (j, k) = exp (2mjk/d) = rf k with r\ = exp {2m /d). (We will make the dependence on d 
explicit below.) In detail 

d-l 

S J , k =J2 F ti>r)A r , k (7) 

r=0 

is a sum of products of scalars times matrices. Since F is invertible, it follows that S is also a basis for the d x d 
matrices. Note that (^) is a special case of (^) with d = 2 and tj = — 1. 

To illustrate these ideas, it is useful to write out the results for d = 3 in detail. Then ?y = exp (27ri/3) and 
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The spin matrices 5 not only form a basis for d x d matrices, but share many other properties with the real Pauli 
matrices, which we record next. We should note that the matrices Sj tk were also defined in an earlier work by Fivel 
Jul on Hamiltonians on discrete spaces, and many of the properties listed below were first established there. 

Proposition 1 Fix d > 2 and let S denote the corresponding set of spin matrices. 

(i) S is an orthogonal basis of unitary matrices with respect to the trace inner product. 

(ii) If d is odd, each matrix in 5 is in SU(d), while if d is even, Sj ik is in SU(d) if and only if j + k is even. 

(Hi) S jtk = (S h oY ■ (S ,i)\ (% fe ) t = V^Sd-j.d-k = {Sj.kY- 1 and [% fc) SV, s ] = {v kr ~ rf s ) S j+r , k+s using addition 
mod d. 

(iv) tr (Sj, k ) = for all (j, k) ^ (0, 0) . 

Proof. The key observation, noted in pj| , is that the matrices are generated by 5i,o and 5*0,1: Sj lk = (Sx t of ■ (Sa.\) k 
with 5*0,1 • 5*1,0 = 775*1,1 . All of the remaining assertions, including orthogonality, follow from those relations and from 
easy computations. A useful consequence of the manipulations is 

(S jjk ) m = v (i - k)Mm - 1)/2 S mj ,mk. (8) 

Unlike the Pauli matrices, these spin matrices need not be Hermitian; for example, when d — 3 only the identity 
matrix is Hermitian. Thus, when computing the coefficients of a density matrix in these bases, as we do next, the 
Hermitian conjugation notation has to be retained. Note that the very last assertion in Corollary 1 corresponds to 
the usual inequality relating the L 2 magnitude of a Fourier transform and the L 2 magnitiude of the original function. 

Corollary 1 (i) The matrix elements of a d x d density matrix p in the different bases are related by (s) = F* ■ (a), 
(ii) s ,o = 1, Sd-j,d-k = V^s^k and z( F ' s )]fl = Pj,j > °- 

(Hi) Ej, fc \ s jM 2 = d H\P]A 2 and \JJ2j, k \ s jA 2 \jH 3 ,k Ift.fcl 2 ^ VVd- 

Proof. We expand an arbitrary density matrix in the two bases: 



d 

j,k j.k 

where aj j. = Tr (a^ k pj gives a,j^ k = Pjj+k, using addition modrf, and Sj jk — Tr ^5*J fe p^. Then from (]7[), Sj tk = 

Sr=o (j, r) a r ,ki which proves (i). Note that we have to include a complex conjugation in the formula which is 
unnecessary in the d — 2 case, since the Hadamard matrix has real entries. The assertions in (ii) follow from the 
definitions and from the fact that p is Hermitian with trace equal to one. Finally the relations in (Hi) follow from 

Tr ((«)* (s)) = Tr ((a) f (F*) (a)) = dTr ((a) t (a)) = dTr (p 2 ), and from Tr (p 2 ) = \\ > £ fe 1/d 2 = 1/d, 

where the X k are the non-negative eigenvalues of the density p. 

Now let N — di x d 2 x • • • x <4 with di > 2 and with the order of multiplication fixed throughout the discussion. We 
use the underlying and fixed tensor product representation of to define the sets of computational and adjusted 
bases E^ and for N x N matrices as 

i®-^ and Afl=Ag! fei ^...0^, 

where j and k correspond to their 6-tuples and the superscripts in parentheses identify the corresponding di . It follows 
that A. ) — E^L U where the addition of the indices is defined by 

j ® k = (ji + k\ mod di, . . . ,jb + kb mod db). (9) 
The corresponding set of spin matrices S*^ is then defined by (5^1) = (A^) or 



N—i 

r=0 



where = F^ ® • • • ® F^ is the usual tensor product of the Fourier transforms F^ which depend on d k - Since 
we will be taking powers of the 77's, we will use subscripts to denote the dependency of 77 on d k : i]k — exp (2m/dk). 
It is easy to show that an equivalent definition of S*^ is given by 



o[N] 



(f^A®) . (10) 
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Linearity again implies that if p^ is a density matrix on the N x N Hilbert space if ^ with 

jn] _x^jn] a [n] _ 1 sr^JN]qim 

P - 1^ a j,k A j,k — 2^ S J,k b j,k > 
3,k j,k 

then 

(>]) = F*^ • (a™) (11) 

and = ■ Thus we have two different representations for a density matrix p^ , and both of them depend on 

the underlying tensor product representation of if^. 

IV. SIGMA VARIATIONS 

A fully separable density matrix can be represented as a convex combination of tensor products of pure states 
or trace one projections, and we need to represent such d x d projections in a systematic fashion in the spin basis. 
(All projections in this paper are trace one projections.) In the Appendix we show how all trace one projections for 
d = 3 can be represented in a form completely analogous to the d — 2 case, but for our immediate purposes we only 
need to characterize a subclass. The motivation is given by writing the particular d = 2 projections h (<tq ± Ofc) as 

Pk{ r ) = SL=o ((~l) r a k) m , where r = or r = 1. Then Pfc(r) is the average of the cyclic subgroup generated by 
(— l) r Ofc, and since ((— l) r Uk) — coj the key property Pfe(r) • Pk(r) — Pk(r) reduces to an exercise in group theory. 
The generalization of this idea to arbitrary d is immediate, and we first treat the case when d is prime. 

Proposition 2 Let d = p > 2 be prime. Let u = (j, k) ^ (0, 0) denote the index of a spin matrix S u , and let r be 
an integer. Then if p > 2, the matrix 

1 p_1 

P u {r) = - ]T (rfS u ) m (12) 

is a projection with unit trace. The assertion is also valid for p = 2 provided iS u is used in lieu of S u throughout 
when u = (1, 1). 

Proof: A matrix P is a pure state or a trace one projection if it is Hermitian, has trace 1, and P 2 = P. First, 
(v r Su) m is proportional to S m j, m k', consequently, it cannot be proportional to So,o for < m < p. Therefore, only S^o 
contributes to the trace of P(u,r), confirming the trace condition. Using (||) it follows that ((ri r S u ) m )^ = (r] r S u ) p m 
and that when p is odd (r] r S u ) m {ri r S u ) p ~ m = {S u f = (rf k ) 2L ^ — S ofi = S ,o. Thus P(u,r) is Hermitian. The 

verification that P (u, r) — P (u, r) follows from an easy computation. The assertion that (rf ) 2 =1 fails for 
prime p only when p = 2 and j = k = 1. Thus, the reintroduction of i and of — cr y = iS^i is required to complete the 
proof. 

As an example of the notation, it is easy to check that if k — 0, then Pj,o(r) is one of the diagonal projections Ei^. 
Other projections are less sparse, however. For example, when d = 3 and k ^ 0, fj,fe(r) has no zero entries in the 
computational basis representation. 

In the preceding proof, we exploited the fact that for d an odd prime the powers of each matrix S u ,u ^ (0,0), 
form a cyclic subgroup of order d. When d is not prime we can get analogous results using a similar proof, but there 
are restrictions on the indices that arise since the coefficient of the identity matrix in (^|) when m = d need not be 
unity. In Proposition 2 this led to the introduction of the factor i = exp (iri/2) when d — 2, and that modification is 
a special case of a more general situation. 

Proposition 3 Suppose d is composite. Let u = (j, k) be (0, 1), (1, 0) or else an index such that j ^ and k ^ 

d-i 

have no common factors. Suppose d is odd or j ■ k is even. Then if r is an integer, P u (r) = h 2 {rf Su) m is a 

m=0 

projection with unit trace. If d is even and j ■ k is odd, then P u {r) = 4 Ylm=o { ar f Su) m is a projection with unit 
trace, where a = e K% l d . 

Proof: Suppose (rf S u ) m or (arf S u ) m is proportional to Sq,o f° r < m < d, so that mj — rd and rak = sd for 
some integers r and s. Since j and k are relatively prime, there are integers a and b such that aj + bk — 1 p0| , and 
it follows that m — ard + bsd = (ar + bs) d, contradicting m < d. Thus P u (r) has trace one. Using (||) when m = d, 
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we find that the coefficient of 5o,o is one in the first case, while in the second case the extra factor of a d = (—1) is 
necessary to make the overall coefficient equal to one. In both cases it follows from that key result as in Proposition 
2 that P^{r) = P u {r) and that P«(r) is Hermitian, completing the proof. 

An important relationship between these subgroup projections and the generating spin matrix follows from the 
definitions. 

Corollary 2 For any integer t > and any d > 2, 

d-l 

(v r S u y=^2v- mt Pu(m + r), (13) 



m=0 



subject to the usual caveat about a. In particular, Sq.q — X)m=o Puijn + r ) 

Proof. £™Jo V - mt P u {m + r) = i E fc (^S u ) k £ m ^-"V* = tfsj, as requ ired. 

Next consider a Hubert space that is the direct product of b Hilbert spaces with dimensions d\, . . . , df,. Projections 
in the constituent di dimensional spaces also define projections in tensor product spaces, and the proof of the following 
is immediate. As before, we let the superscript k denote the dependence on dk- 

Corollary 3 Let N = d\ x e? 2 x • ■ ■ x db and let ff^ = ® a= i H ' da ' be an N dimensional Hilbert space. Let u 
denote a ^-dimensional vector of index pairs Ui — {ji, ki) where < ji, ki < di — 1, and let r = (ri, . . . , r&) where the 
r,; are integers. Then if the Pu^ ( r k) are trace one projections on H^ dk \ 

Pu{r)=(&Pi k J{r k ) 
fe=l 

is a trace one projection on H^ N \ provided aS u is used in place of S u when d is even and u — (J, k) with j ■ k odd. 
Furthermore, if r\ (r) = Ylk=x rfk anc ^ ^ is a non-negative integer, 

(ri (r) S[ N] ) * = E • • • E ® ^^(h + ru), (14) 
i 1= o i b =o fc=i 

and in particular S$ = X) • • • £ ®!Li PvtHh + r k ). 

i 1= o i b =o 

In order to show separability results for Werner densities, we need to identify a special class of fully separable density 
matrices in the tensor product space (d) of n d-dimensional Hilbert spaces, where N = d n . This approach is 
motivated by results in fll6| and is our final variation on the Pauli cr-matrices. 

Proposition 4 Let d > 2 and let «W = (tti, . . . , u n ) and n") = (n, . . . , r„) denote vectors of indices and values 
as defined in the preceding propositions. Then provided aS u is used in place of S u when d is even and u = (j, k) with 
j • k odd, 

\ m=l / 

is a fully separable density matrix on (d). 

Proof. The assertion is true for n = 1 and suppose it holds for n. Let and A n+1 ^ be given index and 

parameter vectors. Since we require only the n'th and (n + l)'st indices in the proof, we leave the other indices fixed 
and implicit and let p(u n ,r n ) denote p (u^ n \r^ n '\ By the induction hypothesis 

^ d-l 

- E POn, r n + s)® Ptti ] ( r »+i - s ) 
is fully separable. Multiplying out and collecting terms produces expressions of the form 

d-l 

M 1 s m ) mi ® • • • ® {rf n s Un ) mi ] ® (?f - +1 s„„ +1 )™ 2 2 E ^ s(mi ~ m2) - 

s=0 

By the same analysis used earlier, terms with mi = m 2 have an overall coefficient of 1 while all other terms have 
coefficient 0, and that completes the proof of the induction step. 
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V. APPLICATIONS 



We now have the tools to prove a sufficient condition for full and partial separability which extend the results in 
ijbSf . This is a general sufficient condition for full or partial separability, and the results in Jl3] and [EjsJ for N = 2™ 
and N = 3™ respectively showing the existence of a neighborhood of the weighted identity in which every density is 
fully separable follow as corollaries. As usual if ^ will denote an N dimensional Hilbert space which can be written 
as a tensor product: = eg) • • • ® , where the H^ dk<l are dfc-dimensional spaces and N — d\ x d 2 x • ■ ■ x db- 

We define D = (di, . . . , d&) and refer to TI^ 1 ! ® • • • ® as the Z? tensor product version of £TM. Since > 

may be represented as a tensor product space in different ways, the kind of separability to be discussed depends on 
the representation. For example, if N = 3™ and H ^ is represented as the tensor product of n three-dimensional 
spaces, we are discussing full separability. If subsets of the trits are taken together and represented in 3 fe -dimensional 
spaces, we are discussing the corresponding partial separability. By virtue of Lemma 1, we know that the fundamental 
mathematics involved doesn't depend on the order in which the tensor products are taken or which trits are grouped 
together. 

In expressing the condition of the theorem, we use the D spin coefficients to introduce an L\ norm on the space 
of N x N densities, and we will refer to that hereafter as the D spin norm and to the related separability as D 
separability. 

Theorem 1 Let H^ N ' denote an iV-dimensional Hilbert space with N = d\ x c?2 x ■ • • X db- Suppose = 
(g) • • • (g) H^ db \ where the ijl dfc l are d^-dimensional Hilbert spaces. If p is a density matrix on H^ N \ then p is 
D = (di, . . . , db) separable provided 



(i,*0#(o,o) 



.[N] 
S 3,k 



< 1, (15) 

where p has the spin representation J^j k s j!k^f!k defined in term of the D tensor product: S^jj — ®^ =1 S^ k ,. It 

follows that in the set of density matrices on there is a neighborhood relative to D of the random state jtSqq 
in which every density matrix is D separable. 

(i) 

Proof: If di is prime or ji and ki are relatively prime, the factor Sj . k . can be written as a weighted sum of projections 

t- (i) 

as in Corollary 2. If di is composite and the indices ji and ki are not relatively prime, then up to a factor of 77^, oj. k . 
can be written as (jlJSu^ for some u — [ji, ki) with ji and ki relatively prime, and thus S^ k . can also be written 
as a weighted sum of projections. Now since p is a density, either S^jj is Hermitian and thus s^j is real, or S^ N k 
appears m a pair s 

f k ] s\ N k ] +s*W (sgl) T . l n the second case we use ( |l4| ) in Corollary 3 and the preceding comments 
to collect the various factors of 77** together and obtain 



+ (CO * = E • • • E ® ^ ('* + ru) {fa^Tf (l) + faaffT, (I) } 

h=o l b =0 k=l 

d\— 1 db — 1 b 

s ffe| E •■■ E (g)^*+r*){exp (0 + exp(-% fc )»7(0} 



l 1= l b =0 k=l 



where 6j t u denotes the phase of /3j,fcsK' and I denotes the b- vector with components The caveat that ctiS u is in the 

projections Pu} (Ik + r k ) in lieu of S u when di is even and u = (ji, ki) with j^ki odd applies throughout the proof and 
will not be explicitly cited. Since on has magnitude 1, only the phase factor will be affected. Using the last assertion 



in Corollary 3, we can write 



S 3,k 



d\ — l djj — 1 b 

s 15 I E • • ■ E (8) & + rfe ) {1 + cos (6 Jtk - arg ( V (I)))} . 

i 1= Z b =0 fc=l 

Since the expression in brackets is non-negative, the right-hand side is a non-negative multiple of a Inseparable density. 
In the case when SjjJ is Hermitian we derive the same expression with the same conclusion. It follows that p can be 
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written as a convex combination of fully separable densities plus the residual term (l — 



(j,fc)#(o,o) 



\ 1 aim 

J N°0,0 



The 



hypothesis guarantees that the coefficient of -hsj^Q is non-negative, and that completes the proof of D-separability. 

As another application of the machinery, we can prove for prime p that the necessary condition s < (l 
is sufficient for full separability of the generalized Werner density matrix (s) = I + s ■ t. We have N = p n , 

I is the identity, r is the projection defined by the state J-!/^) = J2k=o ano - ^ denotes the n-long repeated 
index k . . .k. Given this special structure we find 



1 - s 



p—1 p i 



j=0 k=0 



p-lp-1 

J - * , * \— - ^[N] 
p" p < — ' / —" 3,k 



££^i 

j=a k=o 



where we have used the modular vector addition defined in (|9|). Computing the spin coefficients gives so,o = 1, 
Sj,m — if m is not a k with < k < p, and otherwise 

s ^ = £ F * U,r)-8(r,Ind) 

r ^ 

where Ind = {r : < r < p}. Using the dot product of the index vectors j • r = ^jk^k mod 



Let ine? (p, n) = {j : X^rLo -? r = ^ m °d P}- Then it's easy to check that s ■ ^ = s if and only if j is in /nd (p, n) and 
that there are exactly p n ~ x such indices. All other s ■ ^ equal zero, and we can write WW ( s ) in the spin basis as 



£ 



exp 



-2ni 



p-i 



V* r ) ) a r,k = - U + £ 



exp 



-27TI 



(j • r) 



WW (s) 



1 - s 



°o,o + 6 



p-i 



£ £* 

j£lnd(p,n) k—0 



[N] 



(16) 



Theorem 2 Let p be prime and TV = p n . Then the generalized Werner density matrix W^ N ' (s) is fully separable 
on ffM (p) if and only if s < (l + p n_1 ) _1 . 

Proof. As shown above, necessarily s < (l +p n_1 ) . Checking the preceding derivation, note that 



p-i 



~ £ E ]j - - £ ^ - £ ^ 

j=0 y ]=0 j£lnd(p,n) 

is a sum of fully separable projections. Taking s = (l + p n ~ 1 ) 1 we can write (s) as 

1 



W ln ' (s) 



1 +P 



n-l 



i=0 jelnd(p,n) y \ k=l / 



j£lnd(p,n) 

For each ^ 0, Ind(p, n) is mapped in a one-to-one manner onto itself by j — * fcj where (kj) = kj r mod p. Thus 



1 



1+p 



n-l 



p-1 



;£% + £ i % 5 + £ 



5 



fc=i 



(17) 



i=0 j£lnd(p,n) \ \ 

But since 

(^,0* ® • • • <8> (5,„, 1 ) fe = tfE^jg = S fcj . s 
for j in Ind(p, n), each j-sum in (IT?]) is fully separable by Proposition 4, completing the proof. 
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p 



i+p- 



It follows for the Werner densities that at the extreme value s = (l + p n 1 ) , J2(j fe)^(o o) 
where the coefficients are based on the decomposition D = (p, . . . ,p) . When p = n = 2, that value is 1, showing 
that the global bound of Theorem 1 is attained. However, for larger n and prime p > 2 the condition Wp^ D < 1 
is too strong for that class, and the special structure of the Werner densities allowed a more refined analysis of 
D = (p, . . . , p)-separability. 

It was shown in the qubit case in |l6| that for each n and given e > 0, there exists a D — (2, ... , 2)-inseparable 
density on H^ 2 1 which has || x d < 1 + e - Thus, for each fixed n the sufficient condition of Theorem 1 is the best 
possible for full separability of qubits. We conjecture that the same is true in general: given any separability vector 
D and e > there exists a £>-inseparable density p with ||/^|| x D < 1 + e. 
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VI. APPENDIX 



By emphasizing selected properties of projections for d — 2, we can obtain a representation of all (trace one) 
projections in spin notation for d > 2. We concentrate on d = 3. To motivate the approach, recall from (^) that 
when d — 2, m 3 = (1) po,o + ( — 1) so that this particular spin coordinate is a convex combination of (+1) and 
(— f ), another way of stating the well-known correspondence between 177,3, the coefficient of a z , and the diagonal of p. 
If p is also a projection, then in the computational coordinates, pj^ — bjbk exp (i (ipj — <fk)), so that fixing 7713 fixes 
Po.o = by and pi^i = b\, and only the phase factor 9 — ipo — fi is unspecified. Using the change of basis formula, the 
two remaining spin coefficients of a projection with prescribed 7773 are thus given in terms of the parameter 9 by 

( mi \_f 8 ,i 1 \f b b ie w \ 

\ -ima J \ J \ 1 -1 J V b bie- w J ' 

where < 9 < 2ix. If we let tk denote the value of Sk,i when 9 — 0, we can rewrite the preceding equation as 
Making the obvious definitions, this gives s = M2 (9) t, and we also find that 

The geometry of this result is that if — I < 7713 < 1, then the remaining spin coefficients in the projections associated 
with 7713 can be identified with the range of a one parameter family of invertible mappings {M2 (9)} acting on t and 
are represented by the intersection of the surface of the Bloch sphere with a horizontal plane at height rris . 

The same pattern of results holds for d — 3. Since S2,o = s * 0, the diagonal of a given p is in one-to-one corre- 
spondence with si t o via the equation s^o = po,o (1) + Pi,i (v 2 ) + /°2,2 (v)- That is, S1.0 is a convex combination of the 
vertices of an equilateral triangle in the complex plane and thus uniquely corresponds to the weights of the vertices, 
weights which are the entries of the diagonal of p. For larger values of d, the geometry is more complicated. For 
example if d — 4, the diagonal of a given p corresponds to two spin coefficients: — f < S2.0 < +1 and si.o which is 
restricted to a rectangle in the complex plane with vertices ± (1 + S2,o) /2 ± i (1 — S2.0) /2. In general the diagonal 
of a density matrix p corresponds to d/2 spin coefficients Sjfl, j 7^ 0, when d is even and (d — 1) /2 spin coefficients 
when d is odd. 

Once Sifi is fixed in the d = 3 case, there are three complex parameters remaining to be specified: 30,1,81,1, 
and S2,i, since the other four spin coefficients are forced by the restriction s^-j^-k — if Sj k . If p is a projection, 

S l s j,fc| 2 = IPj.fel 2 = 3, and thus |s2,o| 2 +2 I s m| 2 = 1j tempting one to look for an analogue of the Bloch sphere 
to represent all densities. However, the normalization arising from tr (p 2 ) = 1 is only a necessary condition on the 
parameters, and examples show it's not sufficient. (See also [|l7|.) Instead we follow the d = 2 paradigm and describe 
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trace one projections associated with a fixed si.q. If p is such a projection, then in the computational coordinates 
p = \u) (u\, where \u) denotes a normalized three vector with Uk = bk& lipk and \bk\ 2 — L Fixing si.q fixes the b^s, 
and it follows from Corollary 1 and the structure of p that 







( j 
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j ( 








n 











where 0^ — ifk — fk+i with addition modulo 27r and with the normalization ^2 Ok = mod 27r. Again letting 
denote the value of Sk,i when the 9k 's are chosen to be zero, this time we obtain a two parameter family of projections 
associated with a given value of s^o- Letting s denote the column vector of parameters, t the column vector with 
components tk, and 9 the 3— vector of phase parameters, we have s — M 3 (9) t, where 



1 



3 





i 
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V 
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v 2 ) \ 




I— 1 


1 




I— 1 


>? 




I— 1 


o 
7/^ 





2 



M 3 (0) = - 1 v 2 V • • 1 ?? v 2 = > J(k,e)S , k . (18) 



If + is defined as component- wise addition, then it is easy to check that {M 3 (9)} also defines an Abelian group 
of invertible mappings, 

M 3 (9) ■ M 3 (0) = M 3 (0 + 0) , 

giving the functional equation f (fc, 9) f (j — k, <f) = f (j, 9 + <t>) in analogy with the corresponding result when 
d — 2. We have thus established a correspondence between all trace one projection matrices with given diagonal and 
the range of a two parameter family of mappings acting on t. (We are endebted to Rasmus Hansen for bringing to our 
attention pl[ | , which contains an analysis of the geometry of the convex space of d — 3 densities. The pre-tranform 
characterization of the projections associated with a given diagonal is similar to the results derived here.) 

In the d — 2 case the choices of m z — ±1 produce special cases of projections, and the same is true when d — 3. 
If si.o is one of the extreme points 1, rj, or rj 2 , then two of the bk's equal zero and all of the Sk,i's equal zero. 

It follows that for r = 0, 1, and 2, | S^o + V r Si,o + (v r Si,o) is a trace one projection, and those are the three 

subgroup projections P\ : o(r). A degeneracy which has no analogue in the d — 2 case occurs when si,o lies between 
two extreme points on an edge. Then exactly one of the 6fc's equals zero, and there is a one parameter family of 
projections associated with si.o- The most interesting cases occur when si,o lies in the interior of the equilateral 
triangle. In particular when si.o = 0, the b^s are equal to l/v3, and by choosing the components of 9 appropriately 
from {0, 27r/3, An/3} we find the remaining subgroup projections P u (r). Thus, our entire analysis of separability in 
the d = 3 case uses only the projections associated with the origin and with the vertices of the equilateral triangle. 
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